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Two  experimental  techniques,  electron  swarm  and  electron  beam,  have  been  applied  to  the  problem 
of  electron  attachment  to  POCl3,  with  results  indicating  that  there  is  a  competition  between 
dissociation  of  the  resonant  POClj*  state  and  collisional  stabilization  of  the  parent  anion.  In  the 
electron  beam  experiment  at  zero  electron  energy,  the  fragment  ion  POCF)  is  the  dominant  ion 
product  of  attachment  (96%),  under  single-collision  conditions.  Small  amounts  (~2%  each)  of 
POCI3  and  Cl"  were  observed.  POCI3  and  POClj  ion  products  were  observed  only  at  zero  electron 
energy,  but  higher-energy  resonances  were  recorded  for  POCI',  Cl",  and  C1J  ion  products.  In  the 
electron  swarm  experiment,  which  was  carried  out  in  0.4-7  Torr  of  He  buffer  gas,  the  parent  anion 
branching  ratio  increased  significantly  with  pressure  and  decreased  with  temperature.  The  electron 
attachment  rate  constant  at  297  K  was  measured  to  be  (2.5±0.6)  X  10-7  cm3  s~\  with  ion  products 
POClj  (71%)  and  POClj  (29%)  in  1  Torr  of  He  gas.  The  rate  constant  decreased  as  the  electron 
temperature  was  increased  above  1500  K.  Theory  is  developed  for  (a)  the  unimolecular  dissociation 
of  the  nascent  POClj*  and  (b)  a  stepladder  collisional  stabilization  mechanism  using  the  average 
energy  transferred  per  collision  as  a  parameter.  These  ideas  were  then  used  to  model  the 
experimental  data.  The  modeling  showed  that  Do(Cl-POClJ)  and  EA(POCl3)  must  be  the  same 
within  ±0.03  eV.  ©  2006  American  Institute  of  Physics.  [DOI:  10.1063/1.2176613] 
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I.  INTRODUCTION 

A  number  of  molecules  attach  low-energy  electrons  to 
form  both  the  parent  ion  and  a  dissociative  product.  SF6  is 
the  most  well  studied  of  this  class  of  molecules  due  to  its 
importance  as  an  electron  scavenger  in  electrical  devices.1 
Electron  attachment  to  SF6  is  rapid  and  produces  mainly  SFJ 
at  low  temperatures  and  energies,  but  a  slightly  endothermic 
SFJ  channel  increases  with  electron  energy  and  gas  tempera¬ 
ture.  The  predominance  of  SFJ  production  occurs  even  at 
low  pressure,  indicating  that  a  very  long-lived  excited  SFj* 
species  exists. 

Phosphoryl  chloride  (POCl3)  likewise  attaches  electrons 
rapidly,  with  dissociative  and  nondissociative  channels,  and 
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has  been  the  subject  of  several  recent  studies.2  4  Zero-energy 
electron  attachment  forms  two  main  ionic  species;  the  parent 


anion, 

e~  +  POCI3  — >  POClj,  AH  =  -  1.6  eV,  (1) 

and  POClj  from  a  nearly  thermoneutral  channel 

e~  +  POClj  POClj  +  Cl,  A H  =  0.003  eV .  (2) 

A  third  channel  is  observed  in  very  low  abundance 

e~  +  POCI3  ->  cr  +  POCl2,  A// =0.11  eV.  (3) 

Two  other  channels  are  more  endothermic: 

e~  +  POCI3  —  Cl  J  +  POCI,  A H  =  0.40  eV .  (4) 

<r+POCi3->POcr+ci2,  ah=  1.53  ev.  (5) 


The  reaction  enthalpies  (at  298  K)  given  in  Eqs.  (1 )— (5 )  are 
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from  GAUSSIAN  03  calculations5  reported  in  Ref.  4  and  should 
be  reliable  within  ±0.1  eV.  Electron  attachment  has  been 
found  to  be  almost  as  efficient  as  that  of  SF6  with  little 
variation  in  the  total  rate  constant  with  temperature."  In  con¬ 
trast,  the  ratio  of  the  product  ions  was  found  to  be  very 
dependent  on  temperature;  increasing  temperature  decreased 
the  amount  of  parent  ion  substantially  in  flowing-afterglow 
Langmuir  probe  (FALP)  experiments  performed  at  1  Torr.2 
Increasing  pressure  leads  to  a  substantial  increase  in  the 
amount  of  parent,  so  that  POClj  is  the  dominant  product  at 
pressures  greater  than  100  Torr,  even  at  elevated 
temperatures.3 

Here  we  examine  several  further  aspects  of  the  attach¬ 
ment  to  POCl3.  The  FALP  apparatus  at  the  Air  Force  Re¬ 
search  Laboratory  (AFRL)  has  recently  been  upgraded  to 
allow  for  measurements  as  a  function  of  electron  tempera¬ 
ture.  Measurements  of  rate  constants  were  performed  as  a 
function  of  electron  temperature  at  four  buffer  gas  tempera¬ 
tures.  The  technique  is  essentially  a  clone  of  that  pioneered 
by  Smith  and  Spancl.6  The  results  allow  the  different  effects 
of  internal  and  kinetic  energy  to  be  examined  in  the  same 
experiment,  much  as  is  done  for  ion-molecule  reactions.7 
(See  Christophorou  and  Olthoff  for  descriptions  of  electron 
swarm  experiments,  from  which  similar  data  may  be 
obtained.1)  The  temperature  dependence  of  the  branching  ra¬ 
tios  was  reported  earlier  over  the  temperature  range  of 
296-552  K  at  1  Torr  in  He  gas,2  and  the  pressure  depen¬ 
dence  from  1  to  675  Torr  in  a  N2  buffer  gas,  at  various 
temperatures.3  Here,  we  extend  the  measurements  made  in  a 
helium  buffer  to  cover  a  pressure  range  from  0.4  to  7  Torr  at 
the  same  temperatures  as  earlier,  and  we  have  examined  the 
dependence  of  the  electron  attachment  rate  constant  on  elec¬ 
tron  temperature.  Further,  an  electron  beam  experiment  was 
performed  at  the  Leopold-Franzens  University  to  determine 
branching  ratios  at  essentially  zero  pressure.  The  latter  ex¬ 
periment  also  includes  an  examination  of  higher-energy  reso¬ 
nances,  yielding  dissociative  products.  On  the  theoretical 
side  of  the  problem,  we  have  analyzed  the  competition  be¬ 
tween  the  two  main  channels  for  both  the  new  and  earlier 
data  using  unimolecular  rate  theory.^-11  The  results  of  this 
analysis  may  be  applied  to  other  situations  in  which  the 
POCI3  undergoes  thermal  decomposition.  The  modeling  also 
places  a  narrow  limit  on  the  difference  between  the  electron 
affinity  (EA)  of  POCl3  and  the  bond  dissociation  energy  of 
POCIJ  (into  POCIj+Cl). 

II.  EXPERIMENT 
A.  AFRL  FALP  apparatus 

The  AFRL  FALP  apparatus  has  been  described  in  detail 
previously  and  here  only  details  important  to  the  present  re¬ 
sults  are  given.12  The  FALP  consists  of  a  laminar,  flowing 
electron-ion  plasma  created  by  an  upstream  microwave  dis¬ 
charge  in  He  (~1  Torr).  For  measurements  in  which  the 
electrons  are  equilibrated  with  the  buffer  gas,  that  is,  the 
electron  temperature  Te  is  equal  to  the  gas  temperature  Tg,  a 
small  mount  of  Ar  (—2%)  is  added  to  the  afterglow  to  elimi¬ 
nate  He  metastable  states  and  He2.  The  electron  density  is 
measured  along  the  length  of  the  flow  tube  with  a  movable, 
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cylindrical  Langmuir  probe.  The  plasma  velocity  is  measured 
by  applying  a  small  pulse  to  the  microwave  discharge  and 
timing  the  propagation  of  the  pulse  along  the  flow  tube  axis 
with  the  Langmuir  probe.  In  the  electron-He+-Ar+  plasma, 
the  electron  density  decays  with  distance  due  to  ambipolar 
diffusion. 

In  the  present  work,  POCl3  gas1’  was  added 
(0.2-2  ppmV)  to  the  plasma  at  a  point  downstream  of  the 
discharge  after  the  plasma  was  thermalized,  which  resulted 
in  a  faster  decay  in  the  electron  density  along  the  flow  tube 
axis  downstream  of  the  inlet.  Eventually,  all  electrons  had 
either  undergone  attachment  or  diffused  to  the  walls  of  the 
flow  tube.14  The  data  were  fitted  to  solutions  of  the  one¬ 
dimensional  rate  equations  containing  the  diffusion  and  elec¬ 
tron  attachment  loss  rates.  The  fits,  with  and  without  the 
attaching  gas  present,  yielded  the  electron  attachment  rate 
constant  ka  and  die  ambipolar  diffusion  coefficient  Da.  A 
quadrupole  mass  spectrometer  at  the  terminus  of  the  flow 
tube  was  used  to  identify  the  ion  products  formed  and  (at  low 
mass  resolution)  to  determine  branching  ratios.  Relative  peak 
heights  were  monitored  as  a  function  of  neutral  reactant  con¬ 
centration  in  order  to  distinguish  primary  from  secondary 
product  ions. 

Pressure  dependences  in  the  thermal  branching  ratios 
were  determined  using  the  above  procedure.  Rate  constant 
measurements  were  made  only  at  pressures  from 
0.8  to  1  Torr,  where  Da  is  small  enough  as  to  keep  the  de¬ 
termination  of  ka  reasonably  accurate  (±25%),  and  where 
electron-He  collisions  in  the  plasma  sheath  around  the  Lang¬ 
muir  probe  are  sparse  enough  as  to  not  violate  significantly 
the  Motl-Langmuir  criteria  for  the  probe  operation.  The 
branching  ratios  were  measured  using  very  low  concentra¬ 
tions  of  POCI3  so  that  secondary  ion  chemistry  is  unimpor¬ 
tant.  Since  the  mass  spectra  are  unrelated  to  the  operation  of 
the  Langmuir  probe,  they  can  be  obtained  over  a  larger  pres¬ 
sure  range  than  can  electron  attachment  rate  constants. 

While  He  quenches  electron  energy  rapidly,  Ar  does  not, 
because  of  the  Ramsauer-Townsend  minimum  in  the 
electron-Ar  scattering  cross  section  at  thermal  energies. 
Therefore,  changing  the  buffer  gas  to  Ar  allows  the  study  of 
electron  attachment  at  Te  values  of  thousands  of  Kelvin.  Te 
may  be  varied  by  changing  the  fraction  of  the  Ar  buffer 
flowing  through  the  microwave  discharge  and  by  adding 
some  He  to  the  afterglow.  Measuring  the  electron  tempera¬ 
ture  is  accomplished  by  taking  the  second  derivative  of  the 
current-voltage  characteristic  of  the  Langmuir  probe,  in  the 
electron  retarding  region,  as  described  by  Span  SI.15  Figure  1 
shows  the  second  derivative  of  the  Langmuir  probe  current- 
voltage  characteristic  in  the  electron  retarding  region  for  dif¬ 
ferent  Te.  Electronic  noise,  the  importance  of  patch  fields  on 
the  surface  of  the  probe,  and  the  variation  of  Te  with  distance 
along  the  flow  lube  result  in  an  estimated  20%  uncertainly  in 
Te  in  our  experiment. 

B.  Innsbruck  electron  monochromator 

High-resolution  electron  attachment  to  POCl3  was  stud¬ 
ied  at  the  Leopold-Franzens  University,  utilizing  a  hemi¬ 
spherical  electron  monochromator  in  conjunction  with  a 
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FIG.  1.  Examples  of  fits  to  the  second  derivative  of  the  Langmuir  probe 
current-voltage  characteristic  in  the  electron  retarding  region  to  determine 
the  electron  temperature  Tr.  The  thick  lines  show  the  fitted  region.  The 
curves  have  been  moved  laterally  for  clarity. 

quadrupole  mass  spectrometer.  A  more  detailed  description 
can  be  found  elsewhere.16,17  The  electrons  were  produced  by 
a  heated  filament  and  were  accelerated  with  a  lens  system 
into  a  hemispherical  electrostatic  field  analyzer,  where  the 
electron  energy  resolution  was  set  for  the  present  measure¬ 
ments  to  a  value  of  60-80  meV.  In  the  past,  we  achieved  at 
best  an  electron  energy  resolution  of  35  meV  with  this 
monochromator.  The  present  settings  were  a  compromise 
between  signal  intensity  and  electron  energy  resolution.  Fol¬ 
lowing  the  hemispheres,  the  electrons  were  accelerated  with 
a  second  lens  system  to  the  desired  energy  and  focused  into 
the  collision  chamber,  where  the  electron  beam  intersected 
an  effusive  beam  of  POCl3  at  right  angles.1'1  Anions  formed 
in  the  collisions  were  extracted  by  a  weak  electric  field  (at 
maximum  200meV/cm)  and  were  mass  analyzed  with  a 
high-resolution  quadrupole  mass  spectrometer  (maximum 
mass  range  of  2048  amu).  The  ion  current  was  amplified  by 
a  channeltron-type  secondary  electron  multiplier  operated  in 
a  pulse  counting  mode  and  was  recorded  by  a  computer.  The 
temperature  of  the  gas  beam  was  about  330  K.  The  residual 
gas  pressure  in  the  main  chamber  was  6.8  X  10"8  Torr.  Dur¬ 
ing  the  electron  attachment  measurements  the  pressure  in  the 
chamber  was  1.3  X  10~7  Torr,  measured  with  an  ionization 
gauge. 

The  electron  energy  scale  was  determined  by  measuring 
the  ion  yield  of  a  calibration  gas  under  identical  conditions. 
The  following  attachment  reaction  was  used  for  calibration: 

e_  +  SF6  — >  SFg.  (6) 

The  ion  yield  of  SF^  formed  via  reaction  (6)  exhibits  a  sharp 
peak  at  about  0  eV  formed  via  i-wave  attachment  to  the 
corresponding  neutral  molecules.19  The  electron  energy  reso¬ 
lution  was  determined  as  the  full  width  at  half  maximum  of 
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TABLE  I.  Rate  constants  (±25%)  and  branching  ratios  (±5%  points)  for 
electron  attachment  to  POC1,  using  the  FALP  apparatus  in  1  Torr  He  gas. 
The  small  Cl*  ion  signal  (—1%)  is  not  included  (see  text). 


T(K) 

Rate  constant  (cm’s*1) 

POC1,  fraction 

FOCI;  fraction 

297 

2.5  X 10*7 

29% 

71% 

372 

2.2  X  10*7 

20% 

80% 

457 

2.2  X10*7 

6% 

94% 

552 

2.0  X10*7 

2% 

98% 

the  nominal  0  eV  peak.  The  ion  residence  time  during  ex¬ 
traction  is  estimated  at  100-200  /tts. 

ill.  EXPERIMENTAL  RESULTS 
A.  Electron  attachment  rate  constants 

Both  the  FALP  and  electron  beam  experiments  show  that 
POCl3  attaches  electrons  efficiently  for  near  zero-energy 
electrons.  Electron  attachment  rate  constants  ka  were  mea¬ 
sured  with  the  FALP  apparatus  in  1  Torr  He  gas,  for  tem¬ 
peratures  of  300-552  K,  under  conditions  where  Te=Tg ,  and 
are  listed  in  Table  I  along  with  ion  product  branching  ratios. 
For  reasons  given  below,  the  small  Cl*  ion  product  (~1%) 
was  not  included.  The  measured  ka  are,  on  average,  25% 
greater  than  those  reported  in  our  1998  paper,  which  were 
estimated  accurate  to  within  ±25%.  It  is  not  possible  to  pin¬ 
point  the  source  of  the  difference  except  to  say  that  it  is 
difficult  to  measure  fast  electron  attachment  rate  constants 
with  the  FALP  apparatus:  to  satisfy  first-order  kinetics  re¬ 
quirements,  the  concentration  of  reactant  must  exceed  the 
initial  electron  density  by  a  factor  of  10.  The  lowest  initial 
electron  density  (at  the  reactant  port)  that  is  feasible  in  FALP 
experiments  (for  clear  Langmuir  probe  signal)  is  about 
109  cm-3,  so  that  the  minimum  reactant  concentration  is 
about  10'°  cm"3.  Though  relatively  small  (300  ppbv),  cou- 


Te(K) 

FIG.  2.  Rate  constants  for  electron  attachment  to  POCl3  measured  for  vari¬ 
ous  Tf  and  Tg  in  the  FALP  experiment. 
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FIG.  3.  Resonances  observed  in  electron  attachment  to  POCl3  with  the 
electron  monochromator  experiment,  leading  to  POC1J  and  POClj  ion 
products. 

pling  this  with  the  large  attachment  rate  constant  means  that 
the  reaction  goes  to  completion  over  a  distance  of  only  5  cm 
along  the  flow  tube  axis. 

The  FALP  was  also  operated  with  Ar  gas,  which  allowed 
measurements  to  be  made  at  high  Te  at  different  Tg.  The  ku 
are  plotted  in  Fig.  2  and  are  seen  to  rapidly  decrease  with 
increasing  Te.  The  k„  obtained  at  elevated  Te  are  more  diffi¬ 
cult  to  measure  because  the  electron  densities  tend  to  be 
lower  in  the  Ar  afterglow.  The  ku  at  high  Te  shown  in  Fig.  2 
are  estimated  accurate  to  ±50%,  and  the  electron  temperature 
is  estimated  accurate  to  ±20%.  Within  these  uncertainties, 
data  taken  at  different  gas  temperatures  are  similar. 

B.  Higher-energy  resonances 

Electron  attachment  resonances  were  observed  in  the 
monochromator  experiment  for  POC1J,  POCli,  POC1-,  Cl-, 
and  Clj  ion  products,  as  shown  in  Figs.  3  and  4.  The  first  two 
channels  were  active  only  for  near  zero-energy  electrons. 
The  POCl_  product  appeared  only  in  a  resonance  centered  at 
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FIG.  4.  Resonances  observed  in  electron  attachment  to  POCl3,  leading  to 
POC1",  Cl,,  and  Cl~  ion  products. 

6.4  eV.  The  Cl"  product  was  observed  most  strongly  near 
zero  electron  energy,  followed  by  a  very  weak  resonance  at 
1 .8  eV,  and  a  modest  one  at  6.3  eV.  The  very  weak  Clj 
channel  shows  resonances  at  0.5,  1 .4,  and  6.2  eV.  The  reso- 
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FIG.  5.  Experimental  (symbols)  and  modeled  (curves)  branching  ratios 
[POCI3]/ ([POCIJ] + [POClj]) =57  (S+ D)  in  the  buffer  gas  He.  Modeling  pa¬ 
rameters:  -(A£)/hc=180  cm-1  (dashed  curves)  and  90  cm-1  (solid  curves); 
/rigi,i=0.15. 

nance  at  0.5  eV  is  close  to  the  calculated  endothermicity  for 
the  production  of  C1J,  as  given  in  Eq.  (4).  The  close  vicinity 
of  the  resonances  at  6.4,  6.3,  and  6.2  eV  for  the  three  product 
ions  POCl',  Cl-,  and  C1J,  respectively,  indicates  that  a  com¬ 
mon  transitory  negative  ion  (TNI)  subsequently  dissociates 
into  these  competing  channels.  The  excess  energy  between 
the  position  of  these  resonances  and  the  thermochemical 
thresholds  given  in  Eqs.  (3)-(5)  is  distributed  between  the 
kinetic  energy  of  the  fragments  and  excitation  energy  of  in¬ 
ternal  degrees  of  freedom  of  the  molecular  fragments. 

C.  Stabilization  of  the  parent  anion 

At  very  low  pressures,  where  no  collisions  occur  before 
ions  are  detected,  the  electron  beam  experiment  showed  ion 
products  POCIJ,  Cl-,  and  POCIJ,  with  relative  intensities 
POCIJ/CI'/POCIJ  of  about  47:1.2:1.  That  is,  about  2%  of 
the  parent  anions  formed  in  the  collisions  survives  to  the 
detector,  implying  that  a  few  of  the  nascent  parent  anions 
have  either  a  lifetime  longer  than  the  100-200  /is  flight  time 
or  undergo  stabilization  (via  radiative  or  internal  vibrational 
relaxation).  In  the  FALP  experiment  with  helium  buffer  pres¬ 
sures  in  the  range  of  0.4-7  Toit,  significant  stabilization  of 
the  parent  anion  was  observed  upon  electron  attachment.2 
The  POCIJ  branching  ratio  reached  49%  at  6  Torr  at  room 
temperature.  A  lower  branching  ratio  was  measured  at  el¬ 
evated  temperatures  (e.g.,  10%  at  this  same  pressure  at 
552  K).  Earlier  room  temperature  experiments  in  a  nitrogen 
buffer  had  shown  a  branching  ratio  of  90%  at  675  Torr.'5 
These  results  are  qualitatively  explained  by  resonant  capture 
of  an  electron  into  a  continuum  state  of  POCIJ,  which  then 
undergoes  either  autodetachment,  collisional  stabilization,  or 
decomposition.  The  yields  of  the  parent  anion  in  He  and  N2 
buffer  gases  are  shown  in  Figs.  5  and  6,  along  with  modeling 
results,  of  which  more  shall  be  said  later. 
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FIG.  6.  As  labeled  in  Fig.  5,  for  the  buffer  gas  N2,  with  experimental  data 
from  Ref.  3.  Modeling  parameters:  -(A£)/hc  =  125  cm-1; /rigid =0- 1  (dashed 
curves)  and  0.2  (solid  curves). 

D,  Yield  of  fragment  anions 

Let  us  first  address  the  issue  of  the  low  abundance  of  Cl- 
ions  observed  in  both  experiments.  Calculations  show  this 
channel,  Eq.  (3),  to  be  0.11  eV  endothermic  at  298  K, 
though  the  calculations  carry  an  uncertainty  of  that  same 
order.'1  In  an  earlier  AFRL  publication  on  this  topic,  Cl"  was 
said  to  be  a  direct  product  of  attachment  to  POCl3  because 
the  branching  ratio  increased  from  0.6%  to  5%  over  the  tem¬ 
perature  range  of  296-552  K."  In  the  present  work  with  the 
FALP  apparatus,  with  the  helium  buffer  gas  pressure  varied 
between  0.4  and  7  Torr  and  temperature  varied  between  298 
and  552  K,  the  Cl-  ion  intensity  remained  in  the  neighbor¬ 
hood  of  1  %  of  the  total  ion  intensity  and  did  not  vary  in  any 
systematic  manner  with  pressure  and  temperature.  Thus,  the 
present  conclusion  from  the  FALP  work  is  that  the  Cl-  may 
be  a  result  of  impurities,  possibly  introduced  from  prior  work 
with  chlorine  compounds.  An  experiment  carried  out  at  Mon¬ 
tana  State  University  in  a  N2  buffer  gas  from  1  to  675  Torr 
found  that  the  Cl-  ion  intensity  never  exceeds  1%  of  the  total 
ion  signal.3  However,  in  the  electron  beam  experiments  un¬ 
der  single-collision  conditions,  Cl-  was  observed  at  ~2%, 
and  the  latter  ratio  was  independent  of  pressure  variations  of 
POCI3  in  the  main  chamber  from  1.3  to  6X  10-7  Torr.  The 
positions  of  Cl-  resonances  at  higher  electron  energy  were 
also  independent  of  the  gas  pressure.  The  conclusion  from 
these  and  other  tests  with  the  Innsbruck  electron  monochro¬ 
mator  on  cation  appearance  energies  is  that  the  Cl-  is  a  direct 
product  of  attachment  to  POCl3.  Whether  this  difference  be¬ 
tween  the  two  experiments  is  real  and  related  to  collisional 
effects  or  due  to  impurities  remains  unresolved.  It  is  obvi¬ 
ously  difficult  to  draw  definite  conclusions  about  very  small 
ion  products,  especially  when  the  neutral  is  reactive. 

Under  single-collision  conditions,  the  major  ion  product 
of  attachment  is  POCIJ  observed  in  the  Innsbruck  work  with 
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a  branching  ratio  of  96%.  The  favoring  of  the  POClj  channel 
[Eq.  (2)]  over  the  Cl-  channel  [Eq.  (3)]  is  in  harmony  with 
calculations  which  give  EA(POCl2)  slightly  greater  than  the 
known  EA(C1),  and  consequently,  Dj98(Cl-POClJ)  is  pre¬ 
dicted  to  be  smaller  than  Dj98(POCl2-Cr)  by  0.1  eV.4  Note 
that  these  channels  differ  only  in  the  location  of  the  electron. 
Often  in  dissociative  recombination,  these  are  attributed  to 
differences  in  the  potential  energy  surface  near  the  attach¬ 
ment  region.  However,  the  presence  of  the  parent  anion  re¬ 
veals  that  the  dynamics  are  mainly  controlled  by  unimolecu- 
lar  decomposition  and  not  the  initial  attachment  step.  In  the 
present  FALP  experiments  in  a  He  buffer,  and  in  earlier  work 
in  a  N2  buffer,  smaller  yields  of  POClj  were  observed,  due  to 
collisional  stabilization  of  the  parent  anion.  The  present 
POClj  branching  ratios  are  given  in  Table  I  for  1  Torr  of  He 
gas.  At  other  He  pressures  the  POClj  branching  ratios  may 
be  deduced  from  the  data  in  Fig.  5,  since  POClj  is  the  only 
other  ion  product  in  the  FALP  work  (aside  from  the  ambigu¬ 
ous  Cl"  product). 


IV.  THEORETICAL  TREATMENT 

The  origin  of  the  pressure  dependence  of  the  branching 

between  POOL  and  POClj  was  attributed  to  the 
2—1  J  * 
mechanism  ~ 


POCI3  +  e~—>  POClj*, 

(7) 

POClj* ->  POClj  +  Cl, 

(8) 

POClj*  +  M  — >  POClj  +  M, 

(9) 

where  the  dissociation  step  (8)  involves  a  long-lived  vibra- 
tionally  excited  POClj*  species.  The  present  analysis  pro¬ 
vides  a  quantitative  investigation  of  this  hypothesis.  Denot¬ 
ing  [POClj]  by  S  (stabilization),  [POClj]  by  D 
(dissociation),  the  rate  coefficient  of  process  (8)  by  kdjss,  and 
the  pscudo-first-ordcr  rate  coefficient  of  the  collisional  stabi¬ 
lization  step  (9)  by  kstab[M],  the  branching  ratio  S/(S+D)  for 
POClj  formation  is  symbolically  written  as 

S  _  fcstabM  _  1  (]0) 

S+D  *stab[A/]  + 

^diss  1  ^diss^stabC^] 

In  the  most  rigorous  treatment,  S/{S+D)  would  have  to  be 
determined  by  the  solution  of  a  two-dimensional  master 
equation  in  the  space  of  the  rovibrational  energy  E  and  the 
total  angular  momentum  (quantum  number  J)  of  the  species 
POClj  .  In  this  treatment  kdiss  would  have  to  be  replaced  by 
the  specific  rate  constants  k(E,J)  for  unimolecular  dissocia¬ 
tion  of  POClj  and  kstab  by  the  collisional  stabilization  rate 
coefficients  ycZ,  where  Z  denotes  the  total  rate  constant  for 
energy  transfer  and  yc  is  a  collision  efficiency  resulting  from 
the  solution  of  the  master  equation  for  multistep  rovibra¬ 
tional  energy  transfer.  We  do  not  have  the  required  detailed 
information  on  k{E,J)  and  on  the  state-to-state  collisional 
energy-  and  angular  momentum-transfer  rates,  determining 
yc.  Therefore,  in  the  following  a  simplified  treatment  is  pro¬ 
posed  which  goes  beyond  that  given  in  Ref.  3  and  resolves 
some  of  the  observed  inconsistencies. 


Equation  (10)  indicates  that  S/(S+D)  depends  on  the 
ratio  of  (tdiss/fcstab[^]-  A  comparison  of  modeled  and  experi¬ 
mental  S/(S+D )  thus  can  lead  to  only  a  single  fit  parameter 
which,  to  a  first  approximation,  is  given  by  the  ratio  of  one 
fit  parameter  in  kdiss  and  one  in  kstab.  (As  a  consequence  of 
the  multistep  character  of  the  system  the  fit  parameter  is  only 
approximately  given  by  the  ratio,  see  below.)  In  our  model, 
we  derive  the  fit  parameters  from  the  measured  S/(S+D). 
The  results  then  are  analyzed  and  compared  with  energy 
transfer  and  dissociation  rates  from  related  reaction  systems. 

A.  Specific  rate  constants  for  dissociation 

Because  of  the  lack  of  more  detailed  information  on  the 
transition  state  associated  with  the  dissociation  process  (8), 
£diss  in  Ref.  3  was  estimated  by  simple  Rice-Ramsperger- 
Kassel  (RRK)  theory  in  the  form  kdiss~  ii(E-E0)/E]s~\ 
where  v  is  “the  frequency  of  the  active  oscillator,”  5  is  the 
number  of  vibrational  degrees  of  freedom,  and  E0  is  the  dis¬ 
sociation  threshold  energy.  In  the  present  work  we  adopt  a 
different  point  of  view.  As  the  process  (8)  is  assumed  to 
correspond  to  a  simple  bond  fission  reaction  without  a  re¬ 
verse  barrier,  it  appears  more  appropriate  to  start  with  k(E,J) 
from  phase  space  theory  (PST)  such  as  relevant  for  an  iso¬ 
tropic  potential  energy  surface  for  process  (8).  The  aniso¬ 
tropy  then  is  accounted  for  by  a  rigidity  factor  /rigid  [see  the 
statistical  adiabatic  channel  model/classical  trajectory 
(SACM/CT)  approach  of  Refs.  8-11].  We  therefore  express 
k(E,J)  by 

k(E,J)=fti&id(E,J)k(E,J)PST.  (11) 

Due  to  the  lack  of  information  on  the  potential,  we  simplify 
Eq.  ( 11 )  by  assuming  that  there  is  only  a  small  J  dependence 
of  k{E,J)  (see,  for  instance,  Refs.  9—1 1).  In  a  first  approxi¬ 
mation,  we  further  assume  that  fr\sid(E,J)  is  a  constant  fit 
parameter  over  the  relevant  energy  range.  This  assumption  is 
later  called  into  question.  We  therefore  estimate  kdm  to  be  of 
the  form 

*dta  ~  k(E,J=  0)  ~  /rigidA:(£,i  =  0)PST,  (12) 

where  k(E,J-0)PST  is  obtained  by  convolving  the  number 
of  “open  channels”  (or  “activated  complex  states”)  W(E,J 
_0)psT  for  transitional  modes  of  an  atom + spherical  top 
system  with  the  contributions  from  the  conserved  modes  of 
POClj  (see  Table  II).  We  employ  a  state  counting  algorithm 
such  as  described  in  the  Appendix  of  Ref.  21. 

In  order  to  start  the  fitting  procedure,  a  semiquantitative 
estimate  of/rigid  has  to  be  made.  The  dissociation  process  (8) 
at  long  range  leads  into  a  comparably  weak  ion-induced  di¬ 
pole  potential  which  is  isotropic  and  corresponds  to/rigid=l. 
It  is  known  that  ion-molecule  potentials  of  this  type  at  short 
range  are  governed  by  valence-type  potentials  (see,  for  in¬ 
stance,  Refs.  22-24)  which  correspond  to  values  of  /rigid 
smaller  than  unity.  We  have  analyzed  valence  potentials  from 
ab  initio  calculations  for  neutral  molecules  such  as  H02, 
N02,  and  H202  (see  Refs.  25-27)  and  found  that  /rigid 
(E,J= 0)  was  of  the  order  of  0.1-1  and  showed  only  a  weak 
energy  dependence.  On  the  other  hand,  modeled  short-range 
valence/long-range  ion-induced  dipole  switching  potentials 
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TABLE  II.  Molecular  parameters  used  in  the  modeling. 


Molecule  and 

Frequencies” 

Product  of 

point  group 

(cm-1) 

rotational  constants” 

POCl3(C3„) 

1318.9,  567.4,  567.4,  461.4, 

ABC=2.15  X  10-4  cm-3 

325.7,  325.7,  258.9, 

183.1,  183.1 

POCI3(C,) 

1245.8,  405.6,  387.5,  271.7, 
210.5,  174.8,  131.1,  81.7,  40.1 

POCl;(Ct) 

1214.7,  396.7,  333.3, 
253.6,  201.2,  128.6 

ABC=8.04X  10-4  cm-3 

Additional  quantities:11 

EA(P0C1.,)  =  I.59  eV=hcl2790  cnT1 

£0(POCl;->POCl 
a(He)= 0.205  A3 
a(N2)=  1.740  A3 

;  +  Cl)=hcl2  905  cm-1 

“Frequencies  and  rotational  constants  were  obtained  from  unsealed 
B3LYP/6-31 1+G(34/)  calculations. 

bEA  and  E0  from  G3  calculations  in  Ref.  4;  polarizabilities  from  Ref.  33. 


for  the  dissociations  of  the  cations  of  ethyl-,  n-propyl-,  and 
« -butyl-benzene  in  Refs.  9-11  led  to  energy-dependent /rigid 
consistent  with  experimental  determinations  of  k (£).8  In  this 
case,  /rigicI  was  shown  to  be  of  the  form 

/rigid (E,J)  ~  [1  +  {[E-E0(J)]/Eswrr  (13) 

[with  the  bond  dissociation  energy  E0(J)  and  the  parameters 
£sw/hc=140  cm-1,  m- 2,  and  n-- 2/3  for  C8H|0].10  In  the 
absence  of  more  detailed  knowledge  on  the  potential  of  the 
POCI3  system,  we  started  with  a  constant  /rigid  but  then  also 
tested  an  energy-dependent /rigid  of  the  form  expressed  in  Eq. 
(13). 

Figure  7  shows  k(E,J- 0)  (for/rigid=0.15)  together  with 
thermal  rovibrational  energy  distributions  f(E,T)  of  neutral 
POCI3.  Following  Refs.  2-4,  we  assume  that  during  electron 
attachment  (7)  the  rovibrational  energy  of  POCI3  is  trans¬ 
ferred  into  POCI3  and  is  fully  available  for  dissociation  (8) 
[restrictions  by  centrifugal  barriers  E0{J)  and  additional  con¬ 
tributions  from  the  energy  of  the  electron  are  neglected].  It 
should  be  noted  that  the  rovibrational  average  thermal  energy 
of  POCI3  is  equal  to  <£'th(7’))/hc=990,  1400,  and  1700  cm-1 
for  7=300,  373,  and  423  K  such  as  indicated  by  arrows  in 
Fig.  7. 

B.  Collisional  energy  transfer  rates 

The  overall  rate  coefficients  for  collisional  energy  trans¬ 
fer  between  POCI7*  and  the  buffer  gas  M  are  estimated  to  be 
given  by  the  Langevin  rate  constants  kL=2ve\la/ /x,  with  the 
electronic  charge  e,  the  polarizability  a  of  the  buffer 
gases,  and  the  reduced  mass  /x  of  the  POC1J+M  pairs.  This 
gives  kL(POCl^+ N2)  =  6.4  X  10-10  cm3  molecule-1  s-1  and 
/:ir,(POCl3  +  He)  =  5.4x  10-10  cm3  molecule-1  s_1.  Instead  of 
continuous  state-to-state  collisional  transition  probabilities, 
we  use  a  stepladder  model  with  down  steps  only  where  the 
step  size  is  given  by  the  absolute  value  of  the  average  total 
energy  |{A£')|  transferred  per  collision.  For  highly  excited 
alkyl  benzene  cations,  we  estimated  -(A£)/hc  to  be  of  the 
order  of  200-300  cm-1  for  He  and  N2,  respectively.9,10  As 
(A E)  is  generally  proportional  to  E  and  the  present  energies 
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FIG.  7.  Specific  rate  constants  k(E,J= 0)  (modeled  with  /n£id=0-15)  and 
thermal  distributions  f(E,T)  for  the  dissociation  of  POC1J  (at  T  =  300,  373, 
and  423  K  from  left  to  right  at  small  £;  the  arrows  indicate  average  thermal 
energies  (Elh)  of  POC1J*  if  EA=E0). 

are  about  one-half  of  those  of  the  experiments  with  excited 
alkyl  benzene  cations,9,10  we  estimate  -(A£)/hc  to  be  of  the 
order  of  100-200  cm-1.  We  repeat  that  one  cannot  expect  to 
derive  absolute  values  of  (A£)  from  the  present  measure¬ 
ments  of  S/(S+D).  Instead  one  will  only  be  able  to  derive 
roughly  the  ratios  of  /rigid  and  (A E)  [as  long  as  the  effective 
£diss  is  proportional  to  /rigid  and  kslab  is  proportional  to  (A E) 
in  Eq.  (10),  see  above]. 


C.  Stepladder  model  for  S/(S+D) 

Chemical  activation  systems  such  as  reactions  (7)— (9)  in 
the  past  often  have  been  treated  by  stepladder  models  (see, 
e.g.,  Refs.  28-31)  whenever  more  detailed  information  on 
collisional  transition  probabilities  was  lacking.  Often  solu¬ 
tions  of  master  equations  were  also  found  to  be  equivalent  to 
the  results  from  stepladder  models.31  Applying  a  stepladder 
model  to  the  processes  (7)-(9),  we  have 

S 

S  +  D 


ff  km 

7-7  kL[M]  +  k(E  -  i|<AE)|) 


for  a  given  initial  energy  E.  T(E)  is  equal  to  the  number  of 
steps  between  E  and  E0  and  is  given  by  (£-£0)/|(A£)|  (or 
more  precisely  by  integer  [(£-£0)/|(A£)|]+l).  The  value  of 
S/(S+D)  needs  to  be  thermally  averaged  over  a  thermal 
rovibrational  energy  distribution  f(E,T)  of  POCl3.  We  note 
that  the  energy  E  in  Eq.  (12)  corresponds  to  the  total  rovi¬ 
brational  energy  of  POCIJ  and  that  the  assumption  of  an 
only  weak  J  dependence  of  k(E,J)  only  applies  when  E  is 
defined  in  this  way.  We,  therefore,  finally  calculate 


s  r  T(E) 

—  =  dEf{E,T)H{\+k{E-i\{\E)\)lkL[M]Y\  (15) 

'  +  E>  J  n  1- 1 


where  £=EA+£th  denotes  the  sum  of  the  electron  affinity 
(EA)  of  POCI3  and  its  rovibrational  thermal  energy  £Tth  (E 
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FIG.  8.  As  labeled  in  Fig.  5  with  experimental  points  for  N,  buffer  gas. 
Modeled  curves  are  for  energy  dependent  /ri£id  given  by  Eq.  03)  [full  line, 
E0(J)=Eq(J=  0),  £sw/hc=280  cm"1,  m= 2,  n=-2/3,  -<A£)/hc=125  cm'1], 
for  energy-independent  /rigid  with  EA=£0+0.04  eV  (dashed  line, 
-(A£)/hc=500  cm-1  /ri!,u1=0.15),  and  for  energy-independent  /rigit|  with 
EA=£0-0.04  eV  (dotted  line,  -(A£)/hc=62  cm-1  and  /rigi<)=1.0). 

=0  denotes  the  rovibrational  zero  point  energy  level  of 
POClj). 

D.  Modeling  results 

We  started  our  modeling  by  assuming  that  E0  is  about 
equal  to  EA  (see  Table  II).  Either  varying  (AE),  /rigid  was 
fitted  or  varying  /rigid,  (AE)  was  fitted  such  that  the  experi¬ 
mental  S/(S+D )  curves  as  a  function  of  the  pressure  P  were 
reproduced.  Figures  5  and  6  show  the  results  for  M= N2  and 
He.  One  first  observes  that  the  experimental  temperature  and 
pressure  dependences  within  the  experimental  scatter  are 
well  represented  by  our  model  with  energy-independent  /rigid 
and  temperature-independent  (AE).  Figure  5  for  the  buffer 
gas  M- He  shows  slightly  more  scatter  of  the  measurements 
around  the  modeled  S/(S+D )  curves  than  Fig.  6  for  M=N2. 
Apparently  this  is  due  to  the  larger  uncertainty  of  measuring 
small  POClj  yields.  We  note  again  that  we  are  unable  to 
separate  the  fit  parameters  /rigid  and  (AE)  because  S/(S+D) 
roughly  depends  only  on  the  ratio  /rigid/(A£').  Nevertheless, 
the  fit  of  S/(S+D)  over  wide  temperature  and  pressure 
ranges  appears  quite  satisfactory  and  also  allows  for  extrapo¬ 
lations  outside  the  range  of  the  present  experiments. 

In  addition  to  the  fit  shown  in  Figs.  5  and  6  we  have 
tested  other  possibilities.  First,  we  have  inspected  the  mod¬ 
eling  of  S/(S+D)  with  energy-dependent  rigidity  factors 
such  as  given  by  Eq.  (13),  that  means  with  rigidity  factors 
which  decrease  with  increasing  energy.  We  observed  that  a 
consequence  of  this  increase  in  all  cases  was  a  much  steeper 
increase  of  S/(S+ D)  with  increasing  pressure  than  was  ob¬ 
served  experimentally.  Figure  8  shows  one  example  with 
£sw/hc=280  cm-1,  m= 2,  «=— 2/3  in  Eq.  (13),  and 
-(A£)/hc=  125  cm-1  for  M=N2.  The  experimental  S/(S 


+D)  thus  clearly  point  towards  a  very  weak  energy  depen¬ 
dence  of  /rigjd  such  as  we  observed  for  H02,  N02,  and  H202 
in  contrast  to  the  dissociations  of  cations  of  alkylbenzenes, 
see  above. 

Next  we  investigated  the  dependence  of  the  modeled 
S/(S+D )  on  the  difference  between  the  EA  of  POCl3  and  the 
dissociation  energy  E0  of  POClj.  Figure  8  includes  modeled 
S/(S+D)  for  one  example  with  EA=£'0+0.04  eV  (keeping 
/rigid=0.15  and  fitting  -(A£)/hc=500  cm-')  and  one  ex¬ 
ample  with  EA=£’0-0.04  eV  [choosing  /rigjd=l  and  fitting 
-(AE)/hc=62  cm-1  in  order  to  reproduce  S/(S+D)  at 
300  K],  In  all  cases  either  the  pressure  or  the  temperature 
dependence  of  S/(S+D)  could  not  be  well  represented.  EA 
=E0±0.03  eV  was  about  the  limit  for  which  the  observed 
S/(S+D )  could  be  fitted  over  the  experimental  range.  We, 
therefore,  conclude  that 

EA  =  Eq±  0.03  eV  (16) 

is  about  the  limit  which  is  compatible  with  the  experimental 
S/(S+D).  EA=£0±0.1  eV  was  determined  in  the  G3  calcu¬ 
lations  of  Ref.  4,  in  agreement  with  the  present  evaluations 
of  S/(S+D).  If  EA=  1.60±0.1  eV  is  taken  from  the  G3  cal¬ 
culations  of  Ref.  4  (for  0  K),  being  in  rough  agreement  with 
the  experimental  value  of  EA=1.41  (+0.2/-0.1)  eV  from 
Ref.  32,  then  the  dissociation  energy  of  POClj  in  the  process 
leading  to  POClj+Cl  (for  0  K)  from  the  present  analysis  of 
S/(S+D)  data  follows  as  well  as  £0=1.60(±0.1)  eV. 

We  finally  may  rationalize  our  observation  that  about  2% 
of  the  parent  ions  at  very  low  pressures  have  either  dissocia¬ 
tive  lifetimes  longer  than  100-200  /xs  or  undergo  radiative 
stabilization.  Inspecting  Fig.  7,  one  sees  that  the  low-energy 
tail  of  the  thermal  distribution  reaches  down  to  energies  near 
E0,  where  k(E,J= 0)  approaches  a  threshold  value  of  k(E 
=E0,J=  0)  of  about  100  s-1.  There  is  then  indeed  a  fraction 
of  the  order  of  a  percent  of  the  thermal  distribution  which 
covers  lifetimes  in  the  range  of  100  fx s.  However,  this  re¬ 
quires  that  EA  is  equal  to  £0±0.01  eV,  an  even  more  strin¬ 
gent  restriction  than  postulated  by  Eq.  (16).  One  cannot  rule 
out  that  there  is  also  some  radiative  stabilization,  but  EA 
=£0±0.01  eV  would  be  consistent  with  our  observations. 

A  Cl"  yield  in  the  percent  range  would  also  be  compat¬ 
ible  with  a  pathway  via  POClj.  If  one  shifts  k(E)  by  about 
900  cm-1  to  the  right  in  Fig.  7,  which  corresponds  to  the 
energetics  of  Eq.  (3),  one  has  to  expect  Cl"  yields  from  re¬ 
action  (3)  in  the  percent  range.  However,  our  observations  of 
the  Cl-  yields  are  still  far  too  scattered  and  incomplete  to 
allow  for  more  certain  conclusions. 

V.  CONCLUSIONS 

Experiments  on  the  formation  of  ions  in  electron  attach¬ 
ment  to  POCI3  have  been  described  along  with  a  theoretical 
analysis  of  the  results.  The  experimental  results  consist  of  (a) 
measurements  of  rate  constants  and  product  branching  ratios 
for  electron  attachment  to  POCl3  in  a  He  buffer  gas,  under 
conditions  where  Te-Tg  and  Te >  Tg  (rate  constants  only); 
(b)  determination  of  energies  at  which  resonances  occur  in 
the  attachment  process  for  the  various  ion  products  POClj, 
POClj,  POC1-,  Cl-,  and  Cl  j;  and  (c)  branching  ratios  for  the 
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ion  products  from  single-collision  conditions  up  to  7  Torr  of 
He  gas.  These  final  data  may  be  compared  with  earlier  data 
obtained  at  1-675  Torr  of  N2  gas.1 2 3 4 5  The  theoretical  analysis 
presented  here  treated  (a)  the  unimolecular  dissociation  of 
the  vibrationally  excited  POClj*  formed  in  the  resonant  elec¬ 
tron  attachment  process  with  a  rate  constant  k(E,J)  obtained 
from  phase  space  theory  and  an  anisotropic  contribution  rep¬ 
resented  by  a  rigidity  factor  /rigid  and  (b)  the  stabilization  of 
the  parent  anion  in  collisions  with  the  buffer  gas,  in  terms  of 
a  stepladder  model  using  the  Langevin  rate  constant  coupled 
with  an  average  total  energy  (A E)  transferred  per  collision. 
Modeling  of  the  experimental  data  using  these  concepts  is 
quite  satisfactory,  which  means  that  the  parameters  deter¬ 
mined  for  k(E,J),  /rigid,  and  (A E)  may  be  used  in  other  situ¬ 
ations  where  unimolecular  dissociation  and/or  stabilization 
of  POCI3*  may  be  encountered.  Finally,  the  modeling  of  the 
data  places  a  narrow  limit  of  ±0.03  eV  on  the  difference 
between  D^Cl-POCQ  and  EA(C1). 
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